ADP-01-03/T438 



Non-singlet structure function of the ^He-^H system and 
divergence of the Gottfried integral 

V. Guzey", K. Saito'', M. Strikman^, A.W. Thomas'^ and K. Tsushima" 

■ " Department of Physics and Mathematical Physics, and 

O 

' Special Research Centre for the Subatomic Structure of Matter ( CSSM), 

Q ■ Adelaide University, Adelaide SA 5005 Australia, 

(N 



> 



X 



Tohoku College of Pharmacy, Sendai 981-8558, Japan 
The Pennsylvania State University, University Park, PA 16802, USA 



CN ' Abstract 

o ■ 



We study shadowing and antishadowing corrections to the flavor non- 



^ ■ singlet structure function ^2^*^ ~ -^2^ ^'^'^ show that the difference between 

the one-particle density distributions of ^He and plays an important role 
at very small x. We find that the flavor non-singlet structure function in these 
mirror nuclei is enhanced at small x by nuclear shadowing, which increases 
the nuclear Gottfried integral, integrated from 10~^ to 1, by 11 -i- 36%. When 
integrated from zero, the Gottfried integral is divergent for these mirror nuclei. 
It seems likely that, as a consequence of charge symmetry breaking, this may 
also apply to the proton-neutron system. 
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I. INTRODUCTION 



The measurement of the flavor non-singletQ structure function F2{x,Q^) — -FJ'(x, Q^), 
where Q^) Q^)) is the proton (neutron) structure function, in deep inelastic 

muon-hydrogen and muon-deuterium scattering experiments, performed by the NMC col- 
laboration Ijl]], led to a surprising result. The data revealed an excess of sea down quarks 
as compared to sea up quarks in the free proton. This conclusion has been confirmed by 
the E866/NuSea experiment, where the difference d — u was measured directly using the 
Drell-Yan production of pairs in proton-proton and proton-deuteron collisions 

The results of both experiments contradict the expectation of perturbative QCD (pQCD) 
that u ^ d in the proton. Within the framework of pQCD, the light quark sea is flavor 
symmetric with a good accuracy since it is generated by the perturbative splitting g — > gg, 
which does not distinguish between the u and d flavors. The obvious inconsistency of 
the experimental data with pQCD predictions indicates that non-perturbative effects are 
responsible for creating flavor asymmetry in the light sea quarks. 

The excess of d over u was anticipated well before the measurement on the basis of the 
chiral structure of QCD ||^. Since the NMC experimental discovery and earlier experimental 
indications that d ^ u'ui the proton, this explanation has been actively investigated , with 
the latest discussion centering on the model-independent leading nonanalytic contribution 
IP. Another possible contribution involving the Pauli principle, was flrst explored in pQCD, 
where it was found to give a negligible effect |^. In contrast, non-perturbative calculations 
based on the change in the Dirac sea in the presence of a conflning potential [0 (for recent 
reviews of relevant models also see as well as calculations JTOl, based on the chiral quark- 



^Obviously, the SU(3) flavor symmetry is broken by non-perturbative QCD effects, and the "sin- 
glet" and "non-singlet" combinations of structure functions do not transform as pure singlet and 
non-singlet under SU(3) rotations. In this work, we use the terms "singlet" and "non-singlet" just 
to indicate the quark content of the corresponding structure functions. 
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soliton model of Ref. |jTT|, also predict d > u, which could be of the magnitude observed 
experimentally. Both of these explanations offer considerable insight into the nature of 
hadronic structure in QCD and it is vital to find experimental ways to separate them. 

One way to learn more about the non-perturbative dynamics of the nucleon is to consider 
the non-singlet structure function F^x, Q^) — Q^) for bound nucleons In this 

case, any discrepancy between theoretical predictions and data would indicate that the 
mechanisms, which explained the u ^ d asymmetry for the free proton, are modified in 
nuclear medium. The lightest nuclei which enable one to study the non-singlet combination 



of nuclear structure functions is the pair of mirror nuclei. He and H |12 



The analysis of deep inelastic scattering (DIS) on nuclear targets demonstrates that the 
nuclear environment modifies the properties of the nucleons in a number of ways. At small 
values of Bjorken x, the main effects are nuclear shadowing and antishadowing. In this 
work, we estimate the nuclear shadowing correction to the structure functions, F2^^{x^Q'^) 
for ^He and F^^ix, Q^) for ^H, and for the difference, F2^^{x, Q'^) — F2^{x, Q^), in the region[| 
10~^ < X < 0.02^0.045. The detailed discussion of our approach to the calculation of nuclear 
shadowing is presented in Sect. 0. For larger values of Bjorken x, 0.02 ^ 0.045 < x < 0.2, 
nuclear antishadowing starts to become important. In Sect. we model antishadowing 
by requiring the conservation of the number of valence up and down quarks in '^He and ^H, 
which is a generalization of the baryon number sum rule constraint ||13|| . 

Our results for small x, x < 0.2, can be combined with those of Ref. |jl2[ for the large x 
region in order to present the non-sing let combination {F^^'^ix, Q^)-F2^{x, Q^))/x over the 



full range of Bjorken x. Sect. IV summarizes our results for two models of shadowing and 



^ Note that since the transition region between nuclear shadowing and antishadowing is not 
constrained well by either models or experiments, we use two models of nuclear shadowing with 
different cross-over points between the shadowing and antishadowing regions. This fact is reflected 
in the uncertainty of the upper limit for the shadowing region, x = 0.02 0.045. 
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two pairs of ^He and nuclear wave functions. We also make predictions for the Gottfried 



integral for the ^4=3 system, defined as |T2| 

r3 dx 



/ ^-^{fI^\x,Q')-FI\x,Q'' 



(1) 



In the future, our predictions can be confronted with experiment, for example, with those 
planned at TJLAB Q, RIKEN 111 and RHIC (eRHIC) [111. 



II. NUCLEAR SHADOWING CORRECTION 

The importance of nuclear shadowing in DIS on nuclear targets at small values of Bjorken 
x is experimentally well-established. For recent reviews of the current situation in experiment 



and theory, we refer the reader to Refs. In our approach to nuclear shadowing, we 



choose to work in the target rest frame, where the dynamics of lepton-nucleus interactions 
at small x is transparent. At small Bjorken x, the strong interaction of the virtual photon, 
emitted by the incident lepton, with hadronic (nucleon or nucleus) targets takes place in 
two stages. Firstly, the photon fluctuates into hadronic conflgurations at the distance 
Ic ~ l/{2mNx) before the target 

l7*) = El(^'^l7*)riM , (2) 
it 

where |(/ifc|7*)P is the probability that the photon fluctuates into the state \hk). In pQCD, 
the conflgurations \hk) consist of superpositions of gg, qqg, . . ., Fock states of the virtual 
photon. Secondly, the fluctuations \hk) interact strongly with the target, with some typical 
hadronic cross sections, (7^^^. (Here we have chosen the target to be a nucleus with the 
atomic mass number A.) Within such a picture, the total virtual photon-nucleus cross 
section a^*A can be written: 

'^l*A = Y.\'\^k\l*)?(rhkA ■ (3) 
k 

Here we have suppressed the x and dependence of a^*A for simplicity. The -nucleus 
cross section, (ThkA-, is usually calculated using the high-energy scattering formalism of Gribov 



T8| , which is a generahzation to high energies of the Glauber multiple scattering formalism 



The key element of our approach is an assumption that the sum over the quark-gluon 
fluctuations of the virtual photon in Eq. (|^) can be substituted by some effective state \heff), 
which interacts with bound nucleons of the nuclear target with some effective cross section 
(Teff- Examples of the calculation of nuclear shadowing within such an approximation are 
presented in Refs. [|^,^^-^|. 

In the present work, we will use two models for Ceff- The first model is that of Frankfurt 



and Strikman pSf . The authors used the connection between nuclear shadowing in inclusive 
DIS on nuclei and DIS diffraction in the reaction 7* + p — > X + in order to derive a 
leading- twist model for dg//- Assuming that higher twist contributions to inclusive DIS are 
negligible at Q^=4: GeV^, the model of Ref. ||2^ gives a model-independent]^ description of 
the main contribution to nuclear shadowing (arising from virtual photon scattering off two 
nucleons in the target) in nuclear parton densities and structure functions at small Bjorken 
X. For instance, the leading twist contribution to the nuclear shadowing correction to the 
deuteron structure function F2{x, Q"^) can be calculated unambiguously. For nuclei heavier 
than deuterium, one has to make model-dependent assumptions about cTg// for the scattering 
on three and more nucleons. Since the cross-section fluctuations around the average value 
(Te// practically do not affect shadowing |2^, one can safely use cTg// for the calculation of 
the virtual photon interaction with more than two bound nucleons and employ eikonal or 
quasi-eikonal approximations. (Te// of Ref. p3[ as a function of Bjorken x at Q^=^ GeV^ is 



^ The main assumption which may give a slight model dependence is that the strength of multiple 
rescattering on three or more nucleons is estimated using the quasi-eikonal approximation. Another 
assumption, that the nucleus can be describes as a many-nucleon system, is well justified by 
the small nuclear binding energy per nucleon and also was checked in numerous hadron-nucleus 
scattering experiments at high energies. 
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presented solid line in Fig. 1. 

The second model for cXe//, which we consider, is based on the two-phase model of nuclear 
shadowing for inclusive DIS on nuclei of Ref. [^. This model contains both leading- twist 
(Pomeron and triple Pomeron) and sub- leading twist (vector meson) contributions to Ceff- 
Fig. 1 represents the corresponding (Tg// as a function of Bjorken x at Q^=4 GeV^ as a 
dashed hne. We note that the difference between cTe// of Ref. |^ (solid line in Fig. 1) and 
that extracted from Refs. (dashed line in Fig. 1) lies both in the inclusion of a higher 
twist contribution and in a different parameterization of the Pomeron contribution. 

It is important to note that neither of the models for fXe// distinguishes between virtual 
photon rescatterings on protons and neutrons, i.e. cXe// is a flavor- singlet cross section. In 
this work, we make a simple extension to the fiavor-nonsinglet combination of the virtual 
photon-nucleon cross sections. 

The transition region between nuclear shadowing and antishadowing is poorly known, 
both experimentally and theoretically. In this region, which approximately lies in the range[| 
0.02 < X < 0.07, nuclear structure functions are modified by a host of nuclear effects. Among 
these are nuclear shadowing and antishadowing, two-body nucleon-nucleon correlations in 
the nuclear wave function, the presence of pion degrees of freedom and meson-exchange 
currents. Since our main emphasis is on the very small Bjorken x region, the detailed 
description of the nuclear shadowing-antishadowing transition is unimportant. 

Bearing in mind all these nuclear effects, which if ignored, lead to theoretical uncer- 
tainties in nuclear structure functions, we have included in our analysis the shadowing and 
antishadowing effects only. In addition, we have assumed that the calculations of nuclear 
shadowing, using both models for (Jeff, can be performed most reliably in the range of 
10~^ < X < 0.02. This explains why the upper limit of x in Fig. 1 is set to x=0.02. Since 



^The choice of the lower limit x = 0.02 is motivated by the model of Ref. |2^. The upper limit, 
X = 0.07, corresponds to the largest Bjorken x for which F2'^ /F2 < 1 [25|. 
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(Te// of Ref. p3[ vanishes at x=0.02, we model antishadowing (see Sect. PD in the region 



0.02 < X < 0.2. On the other hand, the two-phase model gives CTg// which is still quite 
significant at a;=0.02 (see Fig. 1). In this case, we force ag// to vanish at a;=0.045 and make 
a linear interpolation between x=0.02 and 0.045. In this case, antishadowing is modelled in 
the region 0.045 <x< 0.2. 

The use of the Gribov-Glauber multiple scattering formalism to calculate CFh^fjA requires 
the |/ie//)-nucleon scattering amplitude and the nuclear wave function. At high energies, the 
|/ie//)-nucleon scattering amplitude fhN{,<i) is purely imaginary with good accuracy. Using 
the optical theorem, fhp{(l) for the proton and fhn{(l) for the neutron are related to the total 



cross sections a^jrj and a^^jj- as 



(4) 



where /3=6 GeV~^ []T^. Here we have assumed that, in general, the effective cross sections 
for the interaction with the proton and neutron are different. 

The ground-state wave functions of ^He and are taken to have a simple Gaussian 



form EOlETl^ 



1=3 



|vl/3HeP oc nexp(-rf /(2«))53(^rl) 
1=1 

|vl/3Hp oc nexp(-rf /(2«'))5=^(E^1) 



(5) 



1=1 



We have checked that the inclusion of the two-body nucleon-nucleon correlations in the 



nuclear wave functions (|g), using the prescription given in Ref. [0], does not change appre- 
ciably the numerical results for nuclear shadowing in the range 10~^ ^ x < 0.05. Hence, we 
shall employ the wave functions of Eq. (|^) in this work. 

The nuclear wave functions of Eq. (^ describe the motion of the centers of the nucleons. 
Thus, the slope parameters a and a' should be chosen to reproduce the nuclear matter radii 
of ^He and '^H. Assuming that only the proton contributes to the nuclear charge radius, the 



nuclear matter radius for a nucleus, R^, takes the form |27] 



Rm - \l Rch - Rl ) (6) 

where Rch and Rp are the charge radii of the nucleus and the proton, Rp = 0.880 ± 0.015 



fm [g8[. In order to estimate the theoretical uncertainty associated with the nuclear wave 
functions, we use two values of the average charge radius of ^He, 1.976 fm and 1.877 fm, 
along with the most recent value of the average charge radius of '^H, 1.76 fm |^9|. From Eq. 
(1), we obtain the following two pairs of matter radii of ^He and '^He: {R^rn"^, -R|^) = (1.769, 
1.524) and (1.658, 1.524) fm. Using the Gaussian-shaped wave functions in the standard 
definition of the average nuclear matter radius, one readily finds that a = This 
leads to the following two pairs of values for the slopes of the nuclear wave functions of ^He 
and (see Eq. (|)): (a, a') = (40.59, 30.06) and (36.11, 30.06) GeV-^. It is important to 
stress that the fact that a 7^ a' is a consequence of the charge symmetry breaking in the 
^He-^H system, which is predominantly the Coulomb repulsion in the ^He system. As will 
be demonstrated later, this leads to the divergence of the corresponding Gottfried integral. 

Using the Gribov-Glauber multiple scattering formalism, along with the elementary scat- 
tering amplitude (^ and the tri-nucleon ground-state wave functions (H), one obtains the 



following I /le//) -nucleus (^He and ^H) total scattering cross sections 



C^3He — ^^eff + ^eff ~t ^ " 



lTi{a + [3) 1447r2(a + [3^ ' 

- a^jf + Za^jf g^^^, ^ e II + ^^^^2^^, ^ . [f) 

Here q\\ = 2mNX is the non-zero longitudinal momentum transferred to the target (with 
the nucleon mass). The negligible x dependence of the triple scattering terms (the last 
terms in the first and second lines of Eqs. (^) is omitted. 

It is convenient to introduce the flavor singlet, a^ff = {cTeff + '^e//)/2) and flavor non- 
singlet, a = a^jrjr — (y^jfi cross sections. Note that the two models of cr^jj and o""^/, those 
of Refs. |2^ and give only the flavor singlet combination [o'^^jf + ^e//)/2- Our analysis 



will demonstrate that the leading contribution of the nuclear shadowing correction to the 
difference -^2'^'^ — -^2"^ is determined by this flavor singlet o"e//- In this new notation, Eqs. 



(0) can be presented as 

C^SHe = ^(^eff + 7:0- — — — 6 H + 



2 87r(a + /5) UAn^ia + 

1_ 3^2^^ - ae//a - 0.25^2 , ^3 - 0.5a2^^a - 0.25ae//a2 + (tVJ 

CTSh = 3cre// cr 7 6 ^11 + 



2 87r(a' + /?) lU7r^{a' + 

It is useful to introduce the short-hand notation 



87T{a + (3) 



1447r2(a + /?)2 ' 



fa' 



-a q 



87r(a' + /3) 



Note that fa, Qa, fa' and (^q,/ are functions of x. Their x dependence originates predominantly 
from the x dependence of CTeff, see Fig. |l]. There is an additional x dependence from the 
non-zero value of gy, which becomes important for x ^ 0.05. 

Using the short-hand notation of Eqs. and ignoring the terms of the order of and 
a^, Eqs. (H) become 

(T3He = 30-e//(l - fa + 9a/^) + |(1 " 2/„ + Qa) , 

CTSH = 3(Je//(l - /a' + fi-a'/S) - ^(1 - 2/a' + S-a') . (10) 

It is important to stress that Eqs. (0) demonstrate that nuclear shadowing in the non- 
vacuum channel (the coefficient in front of the /qO" (/a'O") term) is twice as large as that in 
the vacuum channel (the coefficient in front of the fa<^eff {fa'<^eff) term). This was first 
suggested in Ref. [|l3l. A similar conclusion was reached in the analysis of polarized DIS on 
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He and "^Li [^]. The observation that nuclear shadowing is enhanced by a factor of 



2 in the non-vacuum channel, as compared to the vacuum channel, seems to be a generic 
property of nuclear shadowing and it requires more theoretical work. 
Introducing the structure functions F2{x,Q^) as 



9 



Fi{x, Q') - F^{x, Q2) oc (Tp - a„ = a , (11) 

one can write for the structure functions of ^He and "^H in the shadowing region of Bjorken 

X as 

= 2F^ + - F|(2.5A - g^) - F,^{0.5U) , 

= 2F^ + Fi- F^i2.5U - 9a') - FHO.bU) . (12) 

In Eq. (13), the obvious x and dependence of the structure functions has been suppressed. 
Eqs. ([T^) describe the modification of F2^'^{x,Q'^) and F2^{x,Q'^) at small Bjorken x, as 
a consequence of nuclear shadowing. We observe a qualitatively new effect - the violation 
of SU(2) isospin (charge) symmetry in the wave functions of the A = 3 system, which 
enters through the shadowing correction, induces a violation of SU(2) isospin symmetry for 
the structure functions -^2^'' ^^'^ -^2^- The latter means that ^2^^ and ^2^ are no longer 
related by a rotation in the isospin space. In other words, the charge symmetry violation 
in the wave functions of the A = 3 system results in the SU(2) isospin symmetry breaking 
for nuclear shadowing (regardless of the fact that nuclear shadowing is determined by the 
SU(2)-symmetric exchange with vacuum quantum numbers (Pomeron)). 

As explained above, we assume that using Eq. (p!2D , nuclear shadowing can be calculated 
most reliably in the range 10~^ < a; < 0.02. At higher Bjorken x, nuclear antishadowing 
begins to play a role. Our model-dependent treatment of the antishadowing contribution is 
presented in next section. 

III. NUCLEAR ANTISHADOWING CORRECTION 

The dynamical mechanism of antishadowing is unknown. Thus, at the present stage, 
all considerations of nuclear antishadowing are model-dependent. One possible approach 
to modelling nuclear antishadowing uses the baryon number and momentum sum rules 

10 



T3| , |30|j3T| . The authors of Refs. [|T3| suggest the following scenario, which is consistent with 



the data. (However, it follows from the data only if an assumption is made that higher twist 
effects are small. This assumption is very natural for the case of Drell-Yan data and also 
supported by approximate scaling of DIS data). Nuclear shadowing is present in the valence 
quark, sea quark and gluon parton densities; nuclear antishadowing is present only in the 
valence and gluon parton densities. 

Using Eqs. ([12D we can calculate the nuclear quark parton densities. Adding the anti- 
shadowing contribution to the valence quarks, this leads to: 

= 2d^al + U^al + {2M^al + 0.5M„a/)(-/a + f^d) + dvalQa , 

u^^ = 2u + d+ (2.5m + 0.5J)(-/„) + ug^ , 



/He ^2d + u+ {2.5d + 0.5m) (-/„) + dgo, . (13) 

where u^ai and dyai stand for the valence up and down quark parton densities. The unknown 
functions /^"** and Z^*^** describe nuclear antishadowing for the valence up and down quarks 
in ^He. In order to obtain nuclear quark parton densities in ^H, one needs to replace a by 
a' and m by in the right hand side of Eqs. (pr3|). 

In order to find the functions /^"** and fa"d\ we used conservation of valence up and 
down quarks in ^He 

dxUy^l'ix) = / dx{2u,,ai{x) +4a/(a;)) , 
Jo 

dxd^^i{x) = I dx{2dyaiix) + Uvaiix)) . (14) 
Jo 

The corresponding sum rules are valid for ul^i and dl^i in after the replacement u ^ d 
in the right hand side of Eqs. (0). 

Substituting the first two of Eqs. (|13|) into dH), one obtains the following constraint on 
f^'f and Z,"^"*^ 



a,d 

dx({2.hu.,ai{x) + ^Myai{x))fa - Uyai{x)g^ = f dx(2.5uyai{x) + 0.5dyai{x))f^''^ 

0.0001 ^ Jxo ^ ' 



panti 

5 
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dx(2.5dyaiix) +0.5uyaiix))fa- dyaiix)ga) = [ dx (2.5dyai{x) + 0.5uyaiix)) f^";^\ (15) 

0.0001 ^ ^ JxQ ^ ' ' 



where a;o=0.02 for the calculations with o^jf of Ref. |^ and xo=0.045 for the calculations 
with (Te// based on Refs. [^. In the latter case, we have assumed that a^ff linearly decreases 
from x=0.02 and becomes zero at xo=0.045. This choice of xq is motivated by the NMC 
data on ^He [gl since Ff^jF^ = 1 at x=0.045. 

Using Eqs. ( [TB| ) for ^He and the corresponding equations for '^H, we obtain the following 
equations for the nuclear structure functions 



^He 



2Fl + - F|(2.5/„ - g^) - ^^"(O.S/, 



pp 1 1 1^ ranti ranti\ i rpn / ^ ranti i 57 ranti\\ 

'Q\-'^2val\ g J a,u -^qJ a,d ) ~^ ■'^2val\ e^J a,u -^qJ a,d )J ' 

= 2F^ + F| - F2"(2.5/,, - g^,) - F|(0.5/„0 

1 1 / rpn ( 1 1^ ranti fanti\ , rpp ( ^ ranti , 57 ^(jnti\\ 

where -Ff^^; and F2^^i are the structure functions including only valence quarks. Eqs. (|T6|) 
describe the nuclear shadowing and antishadowing corrections to the nuclear structure func- 
tions Fg^^^ and Fg^^ over the range 10~^ < x < 0.2. 

We would like to stress again that as one can see from Eqs. (|I3|) and (|T^, the violation of 
charge symmetry in the tri-nucleon wave functions induces SU(2) isospin symmetry breaking 
in the quark parton densities and structure functions. In particular, one finds from Eq. (|13D 
that M^^'^ 7^ d^^, and from Eq. ([T6|) that ^2^^*^ is not related to ^2*^ by the permutation 
p ^ n. 

The novelty of Eqs. (|16D consists in the fact that they present the nuclear shadowing 
and antishadowing corrections to structure functions which by themselves are neither fla- 
vor singlet nor flavor non- singlet. Until now, all analyses of nuclear shadowing in DIS on 
nuclei were concerned with nuclei with an equal number of protons and neutrons - i.e., 
flavor singlet nuclei. In applying the previously developed theory of nuclear shadowing and 
antishadowing to DIS on ^He and and deriving Eqs. (0), we have implicitly made the 
following assumptions for the non-singlet combinations of the structure functions F2 and 
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quark densities. We have assumed that Ceff, which controls the amount of nuclear shadow- 
ing is the same for u — and u + d. In other words, in Eq. ([T0|), the same a^ff determines 
the shadowing correction to cxeff (first terms) and a (second term). Another assumption 
was that antishadowing for u and d is the same as for u + d - i.e., it is nil. We believe that 
regardless of the mo del- dependent nature of our estimates, Eqs. ([IB]) provide a reasonable 
estimate of the low x nuclear corrections to the structure functions -^2'^'' and -^2*^. 



IV. RESULTS AND DISCUSSIONS 

Equations (|16D have been used to predict the difference (-Fg^*' — F2^)/x as a function 
of X, in the range of 10~^ < x. In order to test the sensitivity to the input parameters, 
we have considered 5 following combinations of the slopes of ^He and ground-state wave 
functions (||) and models of nuclear shadowing: 

1. a=40.59 GeV-2, a'=30.06 GeV-^, aefj of Ref. [gg with xo=0.02 {xq is the point of 
the transition from shadowing to antishadowing, i.e. fa,u{xo) = fa,d{xo) = /a"**(xo) = 
/m'(^o) = and fa'A^o) = fa'A^o) = fa^^ui^o) = (a^o) = 0. The parameter xq 
enters through Eqs. (|15|)); 

2. a=36.11 GeV-2, a'=30.06 GeV"^, a^fj of Ref. ^ with xo=0.02; 

eff 
eff 



3. a=40.59 GeV'^, a'=30.06 GeV'^, a^ff of Refs. |2§ with a;o=0.045; 

4. a=36.11 GeV-2, a'=30.06 GeV'^, ae// of Refs. H with xo=0.045; 

5. a=40.59 GeV'^, a'=40.59 GeV'^, a^fj of Ref. ^ with xo=0.02. 



For each of these cases, we have assumed the following simple shapes of fa%\ fa4\ fa^,u 
and fa^d- (Below we present only /^'Jj**, with the others being defined in a similar way.) 

f^-f = oml ^^ (^ - ^0)' ^o<x< 0.09 , 

f:f = -^m-x), om<x<xo, 

f^'^^Xx) = 0, elsewhere . (17) 
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The most recent NMC data on F2^'^/F.^ indicates that the antishadowing contribution 
peaks at X = 0.09 |25|. Since DIS on ^He or has not been measured, we assumed that 
antishadowing also peaks at x = 0.09 for DIS on ^He or '^H. This fact is reflected in the 
parameterization of Eqs. (p!7|). The constants h are chosen so that Eqs. (^) are satisfied. 

For quark parton densities in the proton we used the leading order CTEQ5 parameteriza- 
tion CTEQ5L [^. Note also that throughout our work we use the leading order expression 
for the structure functions F2{x,Q'^). This allows us to omit an explicit consideration of 
gluons and forces us to use the leading order quark parton densities, such as for example 
CTEQ5L. 

Figures 2-5 present (^2^^ - F2^)/x as a function of Bjorken x at Q^=4 GeV^ for the 5 
combinations of {a, a') and Ue// given above. The solid lines are results of the calculations 
using Eqs. (O) over the range 10"'^ < a; < 0.2. At larger x, the calculations of Saito et al. 



\ \% were used. For each of the 5 cases, the solid lines should be compared to the dotted lines, 
which present (-Fg^*^ — F2^)/x in the absence of all nuclear effects, when (^2^^*^ — F2^)/x = 
(F| - F^)/x. 

The dotted line in Fig. 2 presents (-^2^^'' — F2^)/x for case 5, when the slopes a and a' 
are chosen to be equal. The large difference between the solid and dotted lines at small x 
demonstrates that the rise of {F2^''-F2^)/x at small X originates from non-cancellation of 
divergent terms in the flavor nonsinglet combination of structure functions F2 of the bound 
proton and neutron, when a 7^ a'. This result implies that charge symmetry breaking (in the 
present case, mainly from the Coulomb force) is very important and enhances the difference 
of the structure functions of mirror nuclei at small x. 

In order to better appreciate the magnitude of these nuclear effects (nuclear shadowing 
and antishadowing) for the flavor non-singlet combinations of structure functions, it should 
be compared to the contribution of nuclear shadowing and antishadowing to the singlet 
combinations of structure functions. Figure 6 presents the flavor singlet combinations F2'^''-|- 
-^2*^ (solid and dashed lines) and 3(^1' -|- F2) (dotted line) as functions of x at Q^=4 GeV^. 
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For the solid and dashed hnes we used the first and third combinations of {a, a') and 
(Te// described in the text aboveQ. One can see from Fig. 6 that, in contrast to the flavor 
non-singlet structure functions, nuclear shadowing decreases -^2'^'^ + -^2*^ as compared to 
3(^2 + F2) but this effect is not as dramatic. The decrease is 4.5% (6%) for the solid 
(dashed) line at x = 10~^. The main conclusion, which one can draw from comparing Figs. 
2, 3, 4, 5 to Fig. 6, is that because of the charge symmetry breaking in the nuclear (^He 
and ^H) wave functions, the nuclear shadowing correction is much more significant for the 
flavor non-singlet combination of structure functions, (-^2'^'' — F2^)/x, than for the flavor 
singlet combination F2^^ + ^2*^. For example, for case 3 and at x = 10~^ and Q^=4 GeV^, 
(^^3He _ p^^nyi^p^Re ^ F2^)=0.0085. However, it should be noticed that since the effect of 
nuclear shadowing in parton densities decreases because of the QCD evolution, we expect 
the ratio (Fg^*^ — ^2^)7(^2^'^ + -^2^) to decrease as increases. 

We used our results for (^2^'^— Fg^)/^ in order to investigated the role played by the small 
X nuclear effects on the Gottfried integral. Table 1 present our estimates of the Gottfried 
integral /^^^''"(lO-^), defined by Eq. (g), and the ratio /5^^''"(10-^)/Jg"(10-4) (/^"(lO-^) 
is the Gottfried integral for the free proton and neutron. We obtained /^.'"(lO^^) = 0.24 
using CTEQ5L, which is in a good agreement with the NMC resuh Jg"(10-^) = 0.235±0.026 
p5| .) We found that the effect of nuclear shadowing increases the Gottfried integral for the 
'^He-'^H system by 11 ^ 36 %, depending on the combination (a, a') and Cg//. 

So far we discussed small, but finite, Bjorken x behaviour of (-^2'^'^ — F2^)/x and the 
integral thereof. At least from the theoretical point of view, one can ask the question: 
what happens to ij^"^' ^(x) when x — >• 0? Our analysis seems to suggest that the Gottfried 
integral for the ^He-'^H system is divergent logarithmically because of the non-cancellation 

^Note also that when o"e// is fixed, the variation of (a, a') leads to very insignificant changes in the 
amount of shadowing and antishadowing. Thus, the solid line in Fig. 6 corresponds to combinations 
1,2 and 5; the dashed line corresponds to combinations 3 and 4. 
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of the factor We observe that this result is not paradoxical since the Gottfried integral 
is not constrained by current algebra - as, for example, the Bjorken sum rule. Thus, the 
value of the Gottfried integral is not related to any physical observable or constant and, in 
principle, can be infinite. 

Our statement that ij^"^' ^(0) diverges is supported by the analysis of the total virtual 
photon-nucleus cross section (structure function F2) at small values of Bjorken x within the 



Gribov model [|18[. Indeed, for DIS on nucleon or nucleus, one can write the dispersion 
integral^ over diffractive masses M for the structure function F2: 

^ ~ 12^ Jo (M2 + g2)2 • (IS) 

Here p(M^) is the the ratio a{e~^e~ hadrons) / a{e^ e~ f^^f^^) {M"^ being mass squared 
of the produced final hadronic state, denoted by ''hadrons"); (t(M^) is the photon-target cross 
section for the production of the final state with mass squared M^. The key assumption 
of the model is that when is constant, x is very small and A ^ 00 (very heavy nuclear 
target), the hadronic fluctuations of the virtual photon \hk) (see Eq. (|])) interact with the 
nucleus with the maximal possible cross section 27iR\ {Ra is the size of the nucleus) [|18[. 
This set of approximations is sometimes called the black body limit. Thus, the contribution 
of the range of for which the black body limit holds, i.e. (j(M^) = 2'kR\, to the structure 
function F2 is 

Q'^2tiR\ i-Mla. dM'^p{M^)M^ 
' ~ 127r3 Jo (M2 + Q2)2 • ( 9) 

The upper limit of integration, M^^x, is defined as the maximal mass squared of a diffrac- 
tively produced intermediate state when, at fixed x and Q^, the black body limit is reached 
for all essential fluctuations of the virtual photon. Within the dipole picture of pQCD, it was 
estimated in Ref. P3| that M^^^ = Q'^xm/x, where Xbu is the critical Bjorken x' entering 



^In general, one has to use the double dispersion representation. However, in the black body limit 
discussed here (see Eq. (|T9|)), only diagonal transitions contribute |18]. 
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the dipole formulation of Ref. [Q, when the black body limit is achieved. The factor Xbbi 



depends on the details of a particular dipole model and, in general, significantly affects the 
absolute value of F2 predicted by Eq. (|1^). Since we are concerned with qualitative and 
model- independent aspects of the x behaviour of F2 following from Eq. (p!9D, after taking 
the integral over the diffractive masses in Eq. (0), we can present the nuclear structure 
function in the form 

F2 oc Q'^ R\\n(l / x) + suhleading terms . (20) 

The application of the black body limit as x — > is also justified for light nuclei and 
nucleons. In particular, using Eq. (]20|), one obtains for the difference of the structure 
functions of ^He and ^H: 

fI^\x, Q^) - F^^{x, Q^) = (3^(i?3He - ^3h) ln(l/a;) + subleadmg terms. (21) 

The charge symmetry breaking in the ^He-^H system manifests itself as the non-equality 
of charge and, hence, nuclear matter radii of ^He (-RsHe) and (-Ran). Substituting Eq. 
(pip into the Gottfried integral yields an integral, divergent as (ln(l/x))^ as x — * 0. Hence, 
we conclude that our analysis of nuclear shadowing and the one within the framework of 
the black body approximation shows that the Gottfried integral for the ^He-^H system is 
divergent. 

It is interesting to note that the above discussed phenomena should also be relevant 
for the free proton and neutron. In this case, similarly to the tri-nucleon system, small 
charge symmetry breaking makes the (hadronic) sizes of the proton and neutron different. 
Specifically, two effects work in the direction of making the radius of the proton larger than 
the radius of the neutron. These are the Coulomb repulsion and the quark mass difference. 
(Since the neutron consists of two d quarks and one u quark and the d quark is heavier 
than the u, the size of the neutron is smaller than that of the proton consisting of two u 
quarks and one d quark.) The difference in sizes should lead to different photoabsorption 
cross sections on the proton and neutron. 
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The analysis of (virtual) photon-hadron interactions demonstrated that in order to suc- 
cessfully describe the data, the photon should contain "soft" and "hard" contributions. The 
soft part interacts with the target with some typical hadronic cross section. Phenomenolog- 
ically, cross sections of soft interactions are proportional to the square of the radius of the 
target hadron - see e.g. Ref. P3. In light of the argument presented above for the size of 
the valence quark distributions in the proton and neutron, the soft component of the photon 
should interact with the proton with a larger cross section. One can expect a similar effect 
for the hard component of the photon. The hard cross section is proportional to the gluon 
field of the target with a cutoff proportional to the size of the target. This makes the cross 
section for interaction with the proton larger than that with the neutron. Hence, in the limit 
of very small values of Bjorken x, the total photoabsorption cross section on the proton is 
larger than on the neutron. In other words, we expect that should be greater than F2, 
which would lead to the divergence of the Gottfried integral, /^'"(O). Further investigations 
of this interesting question are necessary. If, indeed, /^"(O) is infinite, modern parton dis- 
tributions need to be revised since they impose the condition — -FJ* — as x — and, 
hence, give a finite value of Iq"'{0). 



V. CONCLUSION 

We considered the influence of the nuclear effects of shadowing and antishadowing on 
the structure functions -^2'^'^ of ^He and ^2*^ of ^H. We found that these nuclear effects 
increase the Gottfried integral ij^*^' ^(10^'^) by 11 ^ 36%, depending on the model used 
for the nuclear wave functions and for the calculation of nuclear shadowing. We observed 
that the violation of charge symmetry for the nuclear wave functions of ^He and induces 
charge symmetry breaking for the nuclear quark parton densities (as a result of the nuclear 
shadowing correction). This leads to the conclusion that the Gottfried integral, integrated 
over the whole region of Bjorken x, ij^"^' ^(0), is divergent. It is expected that even in the 
case of the free nucleon, the hadronic sizes of the proton and neutron should be different 
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because of the small charge symmetry breaking effect. This suggests that the Gottfried 
integral of the free nucleon should be divergent at very small x. It will be very interesting 
to study the Gottfried integral of the free nucleon at very small x in the future. 

Experiments on DIS off mirror nuclei with large isospin asymmetry should be possible in 
the future [|T3|-|T^. The observation of some deviation from the present calculations would 
provide information on phenomena involving non-pQCD dynamics (like the pion fields) in 
a nuclear medium. If one could vary the atomic number (A) and the difference between the 
proton and neutron numbers {Y = Z — N) independently in measuring the nuclear structure 
functions of unstable mirror nuclei [Q, it would stimulate a great deal of work which may 
eventually lead to genuinely new information on the dynamics of nuclear systems. 
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TABLES 



Case number 




4^"''^(io-^)//g"(io-4) 


1 


0.286 


1.19 


2 


0.267 


1.11 


3 


0.328 


1.36 


4 


0.296 


1.23 


5 


0.242 


1.01 



TABLE L The Gottfried integral Iq^^' ^(10~^), defined by Eq. (|l]), and the ratio of nuclear 
and free space Gottfried integrals Iq^^' ^(10~'^)//^''^(10~^) for the 5 combinations of (a, a') and 
CTe// described in the text. 
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FIGURES 



to 




FIG. 1. Ce// as a function of Bjorken x at Q^=4: GcV'^ from Refs. [23] (solid curve) and [24] 
(dashed curve). Note that the dashed curve includes a higher twist contribution and an earlier 
parameterization of the Pomeron. 
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FIG. 2. Cases 1 and 5. The solid (dashed) hne represents (^2^^ - -^2")/^ ((-^2 - -^2")/^) as a 
function of Bjorken x at (5^=4 GeV^. Case 5 for (-^2^*^ — F2^)/x is given by the dotted hne. For 
parton densities in the proton, the CTEQ5L parameterizations are used. 
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FIG. 3. Case 2. The solid (dashed) hne represents {F2^^-F2^)/x {{F^-F^)/x) as a function 
of Bjorken x at (5^=4 GeV^. . For parton densities in the proton, the CTEQ5L parameterizations 
are used. 



26 



14 



12 



10 




(F2('He).F2('H))/x 
(F,(p)-F,(ii))/x 



10 



-3 



FIG. 4. Case 3. The solid (dashed) hne represents (^2^^ - F2^)/x ((Ff - F^)/x) as a function 
of Bjorken x at (5^=4 GeV^. For parton densities in the proton, the CTEQ5L parameterizations 
are used. 
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FIG. 5. Case 4. The solid (dashed) hne represents {F2^^-F2^)/x {{F^-F^)/x) as a function 
of Bjorken x at (5^=4 GeV^. For parton densities in the proton, the CTEQ5L parameterizations 
are used. 
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FIG. 6. The flavor singlet structure functions F2 + F2 (solid and dashed lines) and 
3(F| + F^) (dotted line) as functions of Bjorken x at (5^=4 GeV^. 
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